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^ ABSTRACT 

^ Time-resolved spectra throughout the orbit of EF Eri during its low accretion 

1^ state were obtained with the Solar Blind Channel on the Advanced Camera for 

^ Surveys onboard the Hubble Space Telescope. The overall spectral distribution 

exhibits peaks at 1500 and 1700A, while the UV light curves display a quasi- 
fy/ sinusoidal modulation over the binary orbit. Models of white dwarfs with a hot 

C/^ spot and cyclotron emission were attempted to fit the spectral variations through- 

^ out the orbit. A non-magnetic white dwarf with a temperature of ~10,000K and 

a hot spot with central temperature of 15,000K generally matches the broad ab- 
sorptions at 1400 and I6OOA with those expected for the quasimolecular H fea- 
^ tures H2 and H^. However, the flux in the core of the Lya absorption does not 

I— I go to zero, implying an additional component, and the flux variations through- 

,__! out the orbit are not well matched at long wavelengths. Alternatively, a 9500K 

^ white dwarf with a 100 MG cyclotron component can flt the lowest (phase 0.0) 

Cf) fluxes, but the highest fluxes (phase 0.5) require an additional source of magnetic 

^ fleld or temperature. The 100 MG fleld required for the UV flt is much higher 

^y~^ than that which flts the optical/IR wavelengths, which would support previous 

O suggestions of a complex fleld structure in polars. 
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1. Introduction 

EF Eri is a well-studied cataclysmic variable that contains a highly magnetic white 
dwarf. It was identified in 1979 as the fourth known Polar (Griffiths et al. 1979; Hiltner et 
al. 1979; Tapia 1979) with an 81 min orbital period. During the next 20 years, it remained in 
a bright high state near 14th mag, with highly modulated optical, IR and X-ray light curves 
due to a high mass transfer rate and X-ray and cyclotron emission from the accretion region 
(Bailey et al. 1982). Models of the cyclotron humps determined magnetic field strengths 
of 16.5 and 21 MG (Ferrario, Bailey & Wickramasinghe 1996). However, in 1997, EF Eri 
entered a low state near 18th mag (Wheatley & Ramsay 1998) which lasted for 9 years, with 
only short excursions to a high state since 2006. Low states are advantageous to study the 
underlying white dwarf and secondary star without the effects of the accretion flow. In the 
low state, the optical spectrum of EF Eri shows broad Balmer absorption lines with Zeeman 
splitting indicating a 14 MG fleld (Wheatley & Ramsay 1998). Beuermann et al. (2000) 
modeled the optical spectrum with a 9500K white dwarf and a 15,000K hot spot over 6% 
of its surface. As there were no features from a secondary star, they concluded that the 
spectral type of the companion must be later than M9. Spitzer mid-IR emission suggested 
the companion was an L or T dwarf (Howell et al. 2006b). Howell et al. (2006a) used the Ha 
emission from this substellar secondary to produce a radial velocity curve, which provided 
a spectroscopic ephemeris as well as a mass estimate for the secondary of 0.055 Mq for a 
white dwarf mass of 0.6 M0. Thorstensen (2003) attempted parallax measurements but the 
faintness of the source and the assumptions used yielded a small range of distances centered 
near 113 pc or a larger range near 163 pc. 

While the accretion rate in low states is reduced by several orders of magnitude from 
the high state, and there is no evidence of a stream of mass transfer during the low state 
of EF Eri, UV studies have shown that hot spots with temperatures of 30,000-70,000K are 
likely present on the white dwarfs of Polars during low states. lUE observations of AM Her 
in the low state revealed an orbital-phase dependent flux modulation, which Heise & Verbunt 
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(1988) interpreted as the signature of moderate heating of a large fraction of the white dwarf 
surface. Gansicke et al. (1995) developed this hypothesis further, showing that a large pole 
cap covering ~ 10% of the white dwarf surface is present both during the high state and the 
low state, and that, at least in the case of AM Her, the luminosity of this heated pole cap 
is comparable to that of the hard X-ray and cyclotron emission from the post-shock region, 
solving the so-called "soft X-ray puzzle" , i.e. the soft X-ray excess with respect to the simple 
reprocessing model of Lamb & Masters (1979) and King & Lasota (1979). Evidence for large, 
"warm" pole-caps was found in many other polars (e.g. Stockman et al. 1994; Gansicke et al. 
2000, Schwope et al. 2002, Araujo-Betancor et al 2005). Gansicke et al (1998) developed a 
3-dimensional model of a white dwarf with a heated pole cap to analyze high-time resolution 
i/ST/GHRS observations of AM Her in a high state, and FUSE and HST/STIS observations 
of AM Her in a low state (Gansicke et al 2006), leading to more precise constraints on the 
size, temperature, and location of the heated pole cap. Especially for the low-state data, 
where no emission from the accretion stream was present, this model provided an excellent 
fit to both the observed modulation of the ultraviolet flux, as well as to the changes in the 
Lyman line profiles. 

GALEX photometric observations of EF Eri during its low state (Szkody et al. 2006; 
2008) in two broad band filters centered near 1550A and 2300A were modeled with a 9500K 
white dwarf and a spot of radius 5.5° with a temperature of 24,000K. Schwope et al. (2007) 
analyzed archival XMM-Newton observations to detect an X-ray flux of 6x10"^^ erg cm~^ 
s~^ during the low state, corresponding to an x-ray luminosity of 2x10^^ erg s~^. While the 
speciflc accretion rate would then be only 0.01 g cm^^ s^^, implying no accretion shock, they 
suggested the X-rays come from some residual accretion at the poles. Li et al. (1994, 1995) 
have produced models for high field white dwarfs where the white dwarf can capture all the 
stellar wind from the secondary and funnel it to the poles. 

Further support for some accretion during the low state comes from the presence of 
cyclotron humps in the IR (Harrison et al. 2004; Campbell et al. 2008). Campbell et 
al. (2008) explored whether high magnetic field cyclotron emission could be producing the 
GALEX UV variations as well as the IR. The high field (240 MG) system AR UMa does 
show UV variability that is likely due to cyclotron emission (Gansicke et al. 2001) and QS 
Tel showed humps at 2200A and 2430A that could be interpreted as cyclotron humps in a 
65 MG field (Rosen et al. 2001). The models of Campbell et al. (2008) showed that the 
UV variations in EF Eri could be cyclotron if the field was 115 MG, while the IR variations 
require a field of 13 MG. Some confirmation of the high field required for the UV is evident 
from the work of Beuermann et al. (2007), who used Zeeman tomography and determined 
that EF Eri has a complex field with some regions as high as 100 MG. 



In order to better understand the cause of the UV hght and variabihty, we obtained 
time-resolved spectroscopy with the Hubble Space Telescope (HST). These data and our 
modehng efforts with a heated polecap and cyclotron emission are described below. 



2. Observations 

The Solar Blind Channel (SBC) on the Advanced Camera for Surveys (ACS) was used 
with prism PRllOL throughout four sequential HST orbits on 2008 January 17. This prism 
provides wavelength coverage with useful sensitivity from ~1200-1900A, the prism results 
in non-linear resolution from ~2A pixel"^ at the blue end to ~40A pixel"^ at the red end. 
The ACCUM mode was used with continuous exposures of 239 sec during each satellite 
orbit. The first orbit had 9 exposures (due to the setup time) while the remaining 3 orbits 
contained 10 exposures. The setup exposure on the target was accomplised with the F140LP 
filter and an integration time of 15 s. The observation times for each orbit are summarized 
in Table 1. 

The reduction package aXel.6 (Kuemmel et al. 2009) provided by STScI was used to 
extract the target and produce a flux and wavlength calibrated spectrum for each 239 sec 
exposure. To obtain the best flux level, we used an extraction width of ±17 pixels cor- 
responding to ± 0.5 arcsec. To show the variability, we created light curves in six bands 
(1220-1260A, 1330-1380A, 1430-1530A, 1530-1640A, 1640-1720A, and 1780-1820A) by in- 
tegrating the fluxes within these regions. Phases were calculated using the spectroscopic 
phasing of Howell et al. (2006a). 

To insure that EF Eri was in its low state during the HST observations, we solicited 
optical ground data from AAVSO members and observatories in the Northern and Southern 
hemisphere in the nights preceding and following the HST times. The New Mexico State 
University (NMSU) Im telescope provided three nights of measuremnts with UBVR filters 
and time-resolved differential light curves in V filter. The Small and Moderate Aperture 
Research Telescope System (SMARTS) obtained B band photometry with a CCD on the 0.9m 
telescope. Spectra throughout an orbit were obtained on the 3.5m telescope at Apache Point 
Observatory (APO) using the Dual Imaging Spectrograph (DIS) with the high resolution 
gratings (resolution ~2A) that provided simultaneous spectra from 3900-5100A and from 
6300A to 7300A. The times of the optical data close to the HST observation are listed in 
Table 1. 



3. Optical Results 

Both the optical photometry and spectroscopy showed that EF Eri remained in its 
deep and extended low state during the HST observations. The calibrated data provided 
B=18.29, V=18.40 and R=18.10 the night preceeding the i/5r observation. The SMARTS 
photometric point from that night is shown in Figure 1 as a star superposed on the B light 
curve accumulated on 68 nights from 2007 June 20 to 2008 March 23 while EF Eri was in a 
low state (Walter 2009). The NMSU V light curves from the nights preceding and following 
the HST times are shown in Figure 2. Both light curves show the normal low state optical 
modulation with a steeper rise to a peak brightness near phase 0.5 and a slower decline to 
minimum near phase 0.0. While the data from January 18th show EF Eri slightly fainter, 
the larger error bars on this night are within the range of variation evident in Figure 1. 
The APO spectra (Figure 3) show the typical low state features of broad Balmer absorption 
lines from the white dwarf flanking weak emission (Wheatley & Ramsay 1998). Howell et al. 
(2006a) determined that the source of the Balmer emission is the substellar secondary and 
is due to stellar activity. The level of emission apparent in Figure 3 is somewhat less than 
that shown in Fig 2 of Howell et al. (2006b) during 2006 January but within the range of 
spectra shown in that paper during 2005. 



4. HST Spectra and Light Curves 

The individual i/ST orbits covered about 0.4-0.5 of the orbital phase of EF Eri on each 
pass, with all phases covered at least once in the 4 orbit sequence. Figure 4 shows the 
variability of the flux within each of the four HST orbits. In this plot, each 239s integrated 
flux point was divided by the average flux, and then converted to a magnitude scale. Due 
to the large variability throughout the binary orbit of EF Eri, we averaged the spectra in 
three ways, which are shown in Figure 5. The solid line is the average of all the data, the 
top dashed line is the average of the 9 integrations within 23% of the peak of the light curves 
(the highest circled points shown in Figure 4), and the bottom dotted line is the average of 
the 9 points within 23% of the minimum flux of the light curves (the lowest circled points in 
Figure 4). The points were selected to avoid the rising and falling portions of the light curve. 
All three averages show the lowest flux at Lya, a broad minimum near 1600A, and peaks 
near 1500 and 1700 A. While these features are preserved in all the averages, the heights of 
the peaks and the depths of the 1600A absorption feature change between peak and trough 
spectra. The 1500A peak becomes less prominent compared to the 1700A peak (which causes 
a corresponding decrease in the depth of the I6OOA absorption feature) during the minimum 
flux times. 



It is possible to interpret these spectra in two ways. If the underlying flux is primarily 
from a cool white dwarf, the I6OOA absorption can be identified as quasi-molecular hydrogen 
H2 which is present in white dwarfs with temperature <13,500K (Koester et al. 1985). There 
is also some evidence for quasi-molecular H^ blueward pf 1400A, which is present in white 
dwarfs cooler than 20,000K. Alternatively, if the UV flux comes primarily from cyclotron, 
the humps at 1500 and 1700A could be caused by cyclotron harmonics. Support for this 
origin is the steep decline longward of ~1750A, which is different from the flatter continuum 
distribution of a cool white dwarf. 

To determine if the drop in flux at longer wavelengths could be a calibration problem, we 
extracted the SBC spectra of two white dwarf standard stars (WD1657+343 and LTT9491) 
taken with the PRllOL and PR130L prisms from the HST archive. The hot white dwarf 
WD1657+343 showed a declining spectrum in both prisms from 1250-1850A while the DZ 
white dwarf LTT9491 showed a discrepancy longward of 1800 A, with the PRllOL spectrum 
being lower than the PR130L spectrum by 5% at 1850 A. To further check, we extracted 
two lUE spectra of LTT9491 from the lUE archive (shown in Figure 6 along with the SBC 
spectrum). This also shows that the SBC calibration is good at short wavelengths but has 
problems longward of 1750A. Thus, we regard the calibrated fluxes of EF Eri to be accurate 
only up to I750A. 

For further phase resolution with good S/N, we created light curves as a function of 
orbital phase by integrating the spectra over the 6 UV bandpasses described in Section 2. 
These light curves are shown in Figure 7. The bandpasses were chosen to span different 
features of the spectra. The 1220-1260A bandpass, which covers the core of Lya, and the 
1780-1820A one, which is the long wavelength downturn in flux, show the least change as a 
function of the orbit. The bandpasses that include the 1500A peak (1430-1530A) and the 
I6OOA absorption feature (1530-1640A) have the largest amplitude variablity (over a factor 
of 2) with peak flux occuring near the same phases as the optical variability (Figures 1 and 
2). However, in contrast to the optical, the UV amplitudes are larger and do not show the 
same asymmetry in shape that is evident in the optical light. 

We attempted to model the spectra and the light curves with hot spots on a white dwarf 
as well as with cyclotron components. 



5. Spot Model 

Using a similar approach to our modeling efforts for AM Her (Gansicke et al. 1998, 
2006), we used the same code to flt the ACS/SBC light curves of EFEri. In brief, the 



surface of the white dwarf is represented by small tiles, with the temperature of each being 
adjustable. The emission of each tile is described by a white dwarf model spectrum with 
the corresponding temperature, and the total emission of the white dwarf is computed by 
integrating the contribution of all tiles on the visible hemisphere. To keep the number of 
free parameters small, the heated region is represented by a circular spot with an opening 
angle ^spot, located at a colatitude /3spot with respect to the rotation axis of the white dwarf, 
and an '^spot with respect to the axis connecting the white dwarf and the donor star. The 
temperature distribution within the spot is assumed to drop linearly in angle from Tccnt at 
the center of the spot to the temperature of the unheated white dwarf, T„d, at the edge of 
the spot. Koenig et al. (2006) modeled cyclotron heating in the case of AM Her, and found 
a spot with a flat temperature near the innermost 10°, which dropped off approximately 
linearly out to ~35°. For modeling the light curves of EF Eri, a linear drop is the simplest 
assumption with the smallest number of parameters. Additional parameters are the white 
dwarf radius -Rwd, and the distance to the system d. We simultaneously fltted the ACS/SBC 
light curves in wavelength ranges 1330-1380 A, 1430-1530 A, 1530-1640 A, and 1640-1720 A. 
The 1220-1260 A and 1720-1820 A light curves were omitted from the flt as they were subject 
to excessive noise and problems in the flux calibration, respectively. Free parameters in the 
flt were Tcs, -Rwd, T'ccnt, ^spot, /Sspot, and ■j/'spot; we kept the distance flxed ai d = 120 pc 
and the inclination of the binary orbit at i = 50°. The free parameters were adjusted 
using an evolution strategy (Rechenberg 1994). It carries out a sequence of mutations of 
the parameters and evaluates them against a fltness function, which we chose to be simply 
X^- The step size of the mutation is itself optimized as a function of the convergence of 
the evolution process. Starting the flt with different initial parameters led to consistent 
convergence in the six-dimensional parameter space. We list in Table 2 the average values 
and standard deviations of six solutions found from very different initial start parameters. 

Inspecting the best-flt light curves to the ACS/SBC observations (Figure 8) reveals that 
the 1330-1380 A data are reasonably well fltted, but that the model signiflcantly over/under 
predicts the observed fluxes in the three longer wavelength bands. This contrasts with the 
quality of the flt achieved for the low-state observations of AM Her (Gansicke et al. 2006). A 
possible explanation is the presence of some residual ultraviolet flux from low-level accretion, 
as suggested by the fact that the flux in the core of Lya (Figure 5) does not drop to zero, as 
expected for the low temperatures of the white dwarf and its pole cap. However, while the 
model fails to exactly match the observed fluxes in the four wavelength ranges, it reproduces 
well the overall spectral shape of the ultraviolet spectra of EF Eri as a function of the orbital 
phase, i.e. as a function of the geometric projection of the heated pole cap (Figure 8). In 
particular, the broad depression seen at all phases near 1600 A corresponds to the quasi- 
molecular H2 absorption that occurs in cool, high-gravity atmospheres (Koester et al. 1985, 



Nelan & Wegner 1985), and the steep rise in flux near 1400 A observed at orbital maximum 
is consistent with the quasi-molecular H^ absorption depressing the flux shortward of 1400A. 
There is some uncertainty as to how the magnetic field affects the lines and quasimolecular 
features in the spectrum of a white dwarf. Gansicke et al. (2001) have shown that the UV 
flux distribution of the highest magnetic field polar AR UMa differs greatly from that of non- 
magnetic white dwarfs in having a flatter UV flux distribution compared to longer optical 
wavelengths. Overall, the low-state ultraviolet spectrum of EF Eri is rather similar to that 
of VVPup (Tcff = 12300 K, Araujo-Betancor et al. 2005), which also exhibits the 1400 A and 
1600 A absorption features. The 1600 A absorption feature disappears in somewhat hotter 
atmospheres, whereas the 1400 A feature remains visible up to ^ 20 000 K (Araujo-Betancor 
et al. 2005, Gansicke et al. 2006). 

The geometry of the heated pole cap is illustrated in Figure 10. The white dwarf 
temperature of ^ 10000 K is consistent with the results from the analyses of optical low- 
state spectroscopy (Beuermann et al. 2000). The white dwarf radius corresponds to a mass 
of ^ O.SMq, but depends on the assumed distance - a true distance larger (lower) than 120 pc 
would imply a correspondingly larger, lighter (smaller, heavier) white dwarf. The maximum 
temperature of the spot, Tccnt — 15000 K is somewhat lower than the value estimated by 
Schwope et al. (2007), which was based on the analysis of the orbital-averaged spectral 
energy distribution of EF Eri as observed with XMM-Newton. Azimuth and co-latitude of 
the heated pole cap are broadly consistent with the location of the main accretion region 
in EFEri as determined by Beuermann et al. (2007) from low-state Zeeman tomography. 
The fractional area of the heated pole cap in EFEri is ~ 22%, which is larger than that in 
AM Her (Gansicke et al. 2006), but comparable to those in V834Cen, BLHyi, MRSer, and 
V895Cen (Araujo-Betancor et al. 2005). 



6. Cyclotron Model 

The cyclotron models were produced using a Constant Lambda (CL) code, which is 
given a more complete overview in Campbell et al (2008). The underlying algorithm is 
based on the assumption that the plasma in the accretion column exists in the large Faraday 
rotation limit >> 1, where = ^^^2^4 Jq neB{s)ds. Radiative transfer therefore decouples 
into two distinct magnetoionic modes: the ordinary (o) and extraordinary (e), whose optical 
depths are given by To^e = ^(t>o,e, where A = lUp/{cUc), where / is the path length through 
the plasma, cOp is the plasma frequency, cOp = {AirNe'^/mY^'^, and coc is the frequency of 
the cyclotron fundamental, Uc = eB/mc, with m being the relativistic mass of the gyrating 
particles. 



Thus, through the parameter A, it is possible to scale a dimensionless emitting slab to 
approximate the observed characteristics of a true 3-D cyclotron emission region. However, 
the model also requires three additional input parameters to calculate the absorption coef- 
ficients (j)o,e'- B, the magnetic field strength, 0, the angle from the magnetic field line to the 
line-of-sight of the observer, and kT, the global isothermal temperature of the emitting slab. 

Due to the large size of the 4-D parameter space, and the computational effort needed 
to enumerate and evaluate all points of the search space, an automated optimization routine 
is required. We used a Genetic Algorithm (GA) to identify the optimal solution (for a good 
primer on genetic algorithms, see Charbonneau (1995)). GAs are a class of optimization 
algorithms inspired by biological evolution, these algorithms encode potential solutions on a 
chromosome- like data structure (called "individuals" ) and apply simulated genetic operators 
to these structures in order to preserve the essential information. GAs begin by establishing 
a random population of individuals, in our case, always respecting the interval of the 4-D 
search space. These individuals then evolve through generations by processes of recombi- 
nation (crossover) and mutation, the natural selection acts at each generation by usually 
prioritizing individuals with higher scores (most well -adapted), given by the fitness func- 
tion. While several publicly available GAs frameworks exist, we chose to use Pyevolve for 
our implementation. Pyevolve is an open-source, modular and object-oriented framework 
for Evolutionary Computation written in Python language by Christian S. Perone (Perone 
2009). Pyevolve is both modular and highly extensible, containing several independent se- 
lection schemes which alter how the fitness function is used to populate the next generation. 

In our work, the individual was constructed out of the four CL parameters with variants 
of each parameter randomly selected from within the region defined by [ 65.00 <B(MG) < 
125.00, 0.1 <kT(keV) < 25, 0.1 <logA < 8, 25° <0< 85°]. A cyclotron model was then 
run for every individual in the population, and then co-added with a white dwarf (WD) at 
120 pc and normalized so that the area under the composite model equals the area under 
an observed spectrum at a specified phase(<l> = 0.00). To enable direct comparison with the 
hotspot results, the WD was computed by using the Gansicke formulation discussed in the 
previous section, but at a lower temperature of 9500 K. The uncertainty in the temperature 
is about lOOOK around a value of 9750K (Schwope et al. 2007). Since the grid of white dwarf 
models in our cyclotron fitting had a spacing of 500K, the 9500K point was chosen as closest 
to the Schwope et al. (2007) value. Our fitness function ranked each individual by looking for 
the minimum the value of x^j normalized to the entire population. We employed a Roulette 
wheel selection scheme, using the Uniform Crossover method with a rate of 90% and a Real 
Gaussian Mutation method with a rate of 15%. It should be noted that convergence to 
global minimum requires maintaining genetic diversity over many generations. To this end, 
a large population of 250 individuals were run against the $ = 0.00 data spectrum over 250 
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generations, finding an optimal solution of xl = 1-24 of B = 100 MG, kT = 6.54 keV, logA 
= 6.88, and = 51.26, which took nearly two weeks of computational resources to evaluate. 

The 39 individual observed spectra of EF Eri were co- added within 0.10 phase bins 
(with about four spectra within each bin) for comparison to the cyclotron models. Orbital 
photometry was estimated by integrating these data through six predefined bandpasses. 
Using an independently determined orbital inclination oi i = 58°, (Campbell et al. 2008), 
and the value of O at cyclotron minimum, O = 51.26°, is sufficient to constrain the magnetic 
co-latitude to /3 = 6.7°. At every other phase-point, a new cyclotron model was computed 
keeping B, kT, and logA constant while allowing G to vary in accordance with the orbital 
geometry determined above. We then integrated the resulting spectra through each band to 
produce our lightcurves. Figure 11 shows the optimized parameter set that was determined 
and the derived quantities while Figure 12 shows cyclotron emission as a function of the 
viewing angle (the cyclotron beam), with the total fiux representing the area under the curve 
shown. The best-fit model at $ = 0.00 is shown in Figure 13 (top). This fit can reproduce 
the minimum spectra well with B=100 MG, kT=6.54 keV, with xt — 1-24. However, the 
corresponding solution at cyclotron maximum produces a rather poor fit to the observed 
spectrum at that phase (Figure 13; bottom). The fit to all 10 phases is shown in Figure 14 
as a stacked series of spectra, and the integrated light is similarly represented as the light 
curves for 6 bandpasses in Figure 15. From all possible fits within a large parameter space, 
we have ruled out a single magnetic component, fixed spot model that can explain the light 
during the entire orbit. Either variable magnetic field strengths or more than one cyclotron 
spot would be needed to account for the variability observed. Given that our use of constant 
A and kT applied to the whole shock is a simplistic representation of a shock that could have 
a variety of temperatures, densities and/or magnetic fields, our work is just a first attempt 
at the problem. 

To further ensure that the Pyevolve was converging on the optimal solution, we also ran 
a large batch of models for B = 100 MG in an evenly spaced grid stepping by (AkT, AlogA) 
= (1.0,1.0) over the ranges of kT= [0:20] and logA = [0:8], respectively. The xl value for 
each were then computed at cyclotron minimum, cyclotron maximum, and also averaged 
over 10 phase-points throughout the orbit. The data were interpolated and mapped to a 
color scheme to produce heat maps shown in Fig. 16. To orbitally modulate the light-curve, 
we again varied to be consistent with an inclination of i = 58°, but computed each set 
of heat maps using various magnetic co-latitudes (/3 = 1, 2, 3 . . . 20°) as inputs. (3 = 7° 
was found to be the optimal solution, where the orbitally averaged heat map had the lowest 
X^. A quick inspection of Fig. 16 demonstrates that the Pyevolve is working. Not only has 
the program converged to the best model at cyclotron minimum (the white spot near kT = 
6.54, logA = 6.88), but an identical solution produces the best model at cyclotron maximum. 
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and indeed averaged through the orbit. Interestingly, while cyclotron minimum has several 
well defined xt minima, the orbital average has no clear best model, presumably because no 
model fits particularly well. 

Finally, we computed the cyclotron fiux, Fcyc, luminosity, Lcyc, and specific mass accre- 
tion rate, m. The following recipe was undertaken: (1) Calculate the accretion spot size. 
Essentially, this is the ratio of the observed cyclotron fiux and the raw model intensity both 
integrated over all wavelengths and subsequently diluted by the surface area of a sphere of 
radius 120 pc. Normalization to the data was accomplished by scaling to lA(Data) which 
is the wavelength integrated model intensity normalized so that the integrated light over 
the range (1220 < Aoata < 1820) equals the integral of the observed fiuxes. lA(Model) is 
just the integrated raw model intensity. (2) The integrated cyclotron fiux was calculated by 
determining the model intensities at all polar angles, and integrating through both angle and 
wavelength. (3) Lcyc is the product of the cyclotron fiux and the spot size. (4) The specific 
mass accretion rate was computed assuming a 0.7 M0 WD, with Kwd determined using the 
Nauenberg (1972) WD mass-radius relationship. The process is also described numerically 
below. 



Wo* = rMD^ta)dA . 

/o A(Model)dAV 120 



00 



Fcyc = 27r / / /A(0)sinedMA (2) 

Jo Jo 

Fcyc^Spot (3) 
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From our final results (Figures 13-15), it is obvious that a single cyclotron component 
can only be fitted to the HST data for a small range in phase (near = 0.00). To fully 
explain the dataset, it is likely that either a second component is necessary, or the addition 
of a hotspot is needed. We are reluctant to explore the implications of a second cyclotron 
component without the inclusion of low-state polarimetry to assess the degree to which the 
UV emission of EF Eri is polarized. 



7. Conclusions 

Our low resolution HST spectra have revealed the first UV spectrum of EF Eri during 
its low accretion state and the corresponding changes in the spectra throughout its orbit. 
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Two interpretations of the lowest orbital UV fluxes at phase 0.0 are possible. If most of the 
flux is from a white dwarf at a distance of 120 pc, the UV spectrum can be matched with 
a temperature of ~10,000K and a hot spot with central temperature of 15,000K (similar 
values found by Beuermann et al. (2000) for the optical region). The broad absorptions 
at 1400 and I6OOA match those expected for the quasimolecular H features H2 and H^ but 
the flux in the core of Lya does not go to zero, implying an additional component. On the 
other hand, the UV spectrum can also be adequately modeled with a 9500K white dwarf at 
120 pc that contributes longward of I6OOA in combination with a high (lOOMG) cyclotron 
component dominating at the shortest wavelengths. 

Both models have difficulty reproducing the observed orbital ffux variations, especially 
for wavelengths longer than 1400A. A normal white dwarf with a hot spot does not produce 
enough amplitude near 1500A but produces too much ffux at 1700 A to account for the 
phase 0.5 modulation. This type of model has problems producing the amplitudes of the 
optical ffux variations as well. The cyclotron models also require some additional component 
(either in temperature or in additional ffeld strengths) to account for the phase 0.5 UV ffuxes 
longward of 1400A. The failure of the simple Constant Lambda model to fft all phases points 
to the need for more complex multi-temperature, multi-lambda models. While our model 
excludes a simple UV emitting pole, a more complex accretion region is very possible. The 
Zeeman tomography work by Beuermann et al (2007) on EF Eri, BL Hyi and CP Tuc shows 
that the ffeld structures are not simple dipoles and there is a range of ffeld strengths in all 
three of these polar systems. For EF Eri, the multipole ffelds range from 4 to 111 MG. 
This range could accomodate the ffelds of 13-21 MG that optimally fft the optical and IR 
spectral regions (Ferrario et al. 1996; Campbell et al. 2008) as well as the UV ffeld strength 
estimated from GALEX (Szkody et al. 2008) and our HST data. An additional uncertainty 
is related to the distance. Ultimately, the choice of model would be best served by UV 
spectropolarimetry to deffne the cyclotron component correctly, but a good parallax and 
UV spectra can pin down the UV luminosity from a high ffeld accreting white dwarf. The 
theoretical proffles for Lya and the quasimolecular H features at high ffelds are not well- 
known, so observations are currently the best way to determine the inffuence of the ffeld on 
the resulting radiation. 

We gratefully acknowledge the help of the AAVSO members who provided magnitude 
estimates close to the time of the HST observations. Support for this research was provided 
by STScI grant HST-GO-11162.01-A. 
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Table 1. Observation Summary 



UT Date Observatory Instrument 



Time 



Int (s) 



2008 Jan 15 


NMSU 


CCD V filter 


02:32-05:34 


300 


2008 Jan 16 


NMSU 


CCD V filter 


01:27-03:09 


300 


2008 Jan 16 


APO 


DIS 


01:48-03:30 


600 


2008 Jan 16 


SMARTS 


ANDICAM B 


02:02:28 


100 


2008 Jan 17 


HST 


SBC 


10:25-11:01 


239 


2998 Jan 17 


HST 


SBC 


11:59-12:39 


239 


2008 Jan 17 


HST 


SBC 


13:35-14:14 


239 


2008 Jan 17 


HST 


SBC 


15:11-15:50 


239 


2008 Jan 18 


NMSU 


CCD V filter 


01:31-04:22 


300 
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Table 2. Best-fit Parameters for White Dwarf plus Hot Spot Model 



Parameter Value 

Teff 9850 ± 175 K 

Tcent 15095 ± 66 K 
i?„d (7.95±8) X lO^cm 
colatitude (/3spot) 17.4 ± 1.4° 

opening angle (6'spot) 57.2 ± 2.2° 

azimuth (V'spot) 4.1 ± 0.4° 
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Phase (P = 0.056 days) 



2.0 



Fig. 1. — SMARTS composite B data for the extended low state during which the HST 
observations were made. The hght curve is comprised of 152 points on 68 nights from 2007 
June 20 to 2008 Mar 23. The star near phase 0.7 is the night of 2008 Jan 16. 
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Fig. 2. — NMSU V data relative to a comparison star for Jan 15 (dots), Jan 16 (triangles) 
and Jan 18 (squares). Error bars are ±0.03-0.05 mag for Jan 15 and 16 and ±0.06-0.08 mag 
on Jan 18. 
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Fig. 3. — DIS blue and red spectra averaged over the orbit of EF Eri. The Balmer emission 
shows that some level of activity on the secondary star is present. 
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Fig. 4. — Variability of EF Eri throughout each i7ST orbit. Each point is a 239s integration, 
with a differential magnitude computed by dividing the integrated flux by the average flux 
over all 4 orbits and converting to a magnitude scale. The 9 circled points near minimum 
light are those used for the minimum flux values (within 23%), and the 9 circled points near 
maximum light are those within 23% of the maximum flux values. 
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Fig. 5. — Average spectrum with an extraction width of 17 using all 39 data points (middle 
solid line), with the 9 points within 23% of the lowest flux points (bottom dotted line) and 
with the 9 points within 23% of the highest flux points (top dashed line). 
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Fig. 6. — Two lUE spectra of LTT9491 (high resolution lines) compared to the SBC spectrum 
(low resolution line). Continuum fluxes match well until 1750A where the SBC flux declines 
relative to the lUE values by about 
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Fig. 7. — Light curves in 6 UV bands as a function of spectroscopic orbital phase. Each 
point is the integration of the flux within the bandpass. 
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Fig. 8. — Best fit to the ACS/SBC light curves in four wavelength bands using the white 
dwarf plus pole cap model. The observed ACS/SBC fiuxes are indicated by different symbols, 
the models corresponding to these bands by different line styles. 
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Fig. 9. — Observed (gray lines) and model spectra at orbital minimuni (bottom curves) and 
orbital maximum (top curves). The model spectra are generated from the best-fit to the 
ACS/SBC light curves (Figure 7). The broad depressions near 1400 A and 1600 A are the 
quasi-molecular H^ and H2 absorptions typical of cool, high-gravity hydrogen atmospheres. 
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Fig. 10. — Geometry of the heated pole cap in EFEri. The unheated white dwarf is repre- 
sented by grid hnes, the pole cap is gray-shaded. Orbital phases are 0.1, 0.3, 0.5, 0.7, and 
0.9 (from left to right). 
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Fig. 11. — Derived orbital geometry, accretion spot size and accretion rates and the orbital 
variation of the viewing angle of the magnetic pole. is the current orbital phase and (po is 
maximum light. 
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Fig. 12. — Dependence of the total integrated intensity as a function of the viewing angle 
G. The area under the curve represents the total cyclotron flux. 
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Fig. 13. — Best fit model spectrum to the phase 0.0 iKT spectrum (top) and phase 0.5 EST 
spectrum (bottom). Data are points with error bars. Lower curves are a 9500K white dwarf, 
mid curves are the cyclotron component for B=100 MG and top curves are the sum of the 
two components. 
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HST SBC Phase-Resolved Spectra of EF Eri 
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Fig. 14. — Phase-resolved HST spectra (black with error bars) along with the cyclotron 
(red), white dwarf (blue) and co-added CYC+WD model (green) for the phase = 0.0 data. 
Each subsequent phase shows only the data and co-added models. Fluxes at phase 0.0 are 
real, each subsequent phase has a constant offset of lxlO~^^ added to the flux for clarity of 
presentation. 
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Fig. 15. — Light curves in 6 bandpasses as a function of orbital phase, with the white dwarf, 
cyclotron and sum model shown along with the 10 phase points created in each bandpass by 
integrating the fluxes within each bandpass region listed on the right. 
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Fig. 16. — x^ heat maps of the kT-logA parameter space for B=100 MG and P=7° at phases 
0.0, 0.5 and the average over the entire orbit. 



